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The behaviour of magnetic field in anisotropic Bianchi type I cosmological model for bulk viscous 
distribution is investigated. The distribution consists of an electrically neutral viscous fluid with an 
infinite electrical conductivity. It is assumed that the component u\ of shear tensor cr^ is proportional 
to expansion (9) and the coefficient of bulk viscosity is assumed to be a power function of mass 
density. Some physical and geometrical aspects of the models are also discussed in presence and also 
in absence of the magnetic field. 
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1. Introduction 

A resurgence of interest in anisotropic, general-relativistic cosmological models of the 
universe 1,2 has been stimulated by the discovery of the cosmic microwave radiation, the 
study of its isotropy, the problem of the abundance of primordial helium, and the possi- 
bility of large-scale primordial magnetic fields. The choice of anisotropic cosmological 
models in Einstein system of field equations leads to the cosmological models more gen- 
eral than Robertson- Walker model. 3 The occurrence of magnetic fields on galactic scale is 
well-established fact today, and their importance for a variety of astrophysical phenomena 
is generally acknowledged as pointed out by Zeldovich et al. 4 Also Harrison 5 has sug- 
gested that magnetic field could have a cosmological origin. As a natural consequences 
we should include magnetic fields in the energy-momentum tensor of the early universe. 
The presence of primordial magnetic fields in the early stages of the evolution of the uni- 
verse has been discussed by several authors. 6 ~ 15 Strong magnetic fields can be created 
due to adiabatic compression in clusters of galaxies. Large-scale magnetic fields give rise 
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to anisotropics in the universe. The anisotropic pressure created by the magnetic fields 
dominates the evolution of the shear anisotropy and it decays slower than if the pressure 
was isotropic. 16 ' 17 Such fields can be generated at the end of an inflationary epoch. 18-22 
Anisotropic magnetic field models have significant contribution in the evolution of galaxies 
and stellar objects. Bali and Ali 23 obtained a magnetized cylindrically symmetric universe 
with an electrically neutral perfect fluid as the source of matter. Several authors 24 ~ 28 have 
investigated Bianchi type I cosmological models with a magnetic field in different context. 



Most cosmological models assume that the matter in the universe can be described by 
'dust' (a pressure-less distribution) or at best a perfect fluid. Nevertheless, there is good 
reason to believe that - at least at the early stages of the universe - viscous effects do play 
a role. 27-29 For example, the existence of the bulk viscosity is equivalent to slow process 
of restoring equilibrium states. 30 The observed physical phenomena such as the large en- 
tropy per baryon and remarkable degree of isotropy of the cosmic microwave background 
radiation suggest analysis of dissipative effects in cosmology. Bulk viscosity is associated 
with the GUT phase transition and string creation. Thus, we should consider the presence 
of a material distribution other than a perfect fluid to have realistic cosmological models 
(see Gr0n 31 for a review on cosmological models with bulk viscosity). The effect of bulk 
viscosity on the cosmological evolution has been investigated by a number of authors in 
the framework of general theory of relativity. 32-46 



In recent years, models with relic cosmological constant A have drawn considerable 
attention among researchers for various aspects such as the age problem, classical tests, 
observational constraints on A, structure formation and gravitational lenses have been dis- 
cussed in the literature. It is remarkable here that in the absence of any interaction with 
matter or radiation this would force the cosmological constant to be constant, but, in the 
presence of the interaction with matter or radiation, a solution of Einstein's field equation 
and assumed equation of covariant conservation of energy with a time varying A can be 
found. For these solutions, conservation of energy requires that any decrease in the energy 
density of the vacuum component be compensated for by a corresponding increase in the 
energy density of matter or radiation. Some of the recent discussions on the cosmolog- 
ical constant "problem" and consequence on cosmology with a time-varying cosmologi- 
cal constant are contained in Bertolami, 47 Ratra and Peebles, 48 Dolgov, 49-51 Sahni and 
Starobinsky, 52 Padmanabhan, 53 Vishwakarma. 54 Recent observations by Perlmutter et al. 
55 and Riess et al. 56 strongly favour a significant and positive A. Their finding arise from 
the study of more than 50 type la supernovae with redshifts in the range 0.10 < z < 0.83 
and suggest Friedmann models with negative pressure matter such as a cosmological con- 
stant, domain walls or cosmic strings (Vilenkin, 57 Garnavichef al. 5S ) Recently, Carmeli 
and Kuzmenko 59 have shown that the cosmological relativity theory (Behar and Carmeli 60 ) 
predicts the value A = 1.934 x 10 -35 s -2 for the cosmological constant. This value of A is 
in excellent agreement with the measurements recently obtained by the High-Z Supernova 
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Team and Supernova Cosmological Project (Garnavich et al.; 58 Perlmutter et al.; 55 Riess 
et al.; 56 Schmidt et al. 61 ) The main conclusion of these works is that the expansion of the 
universe is accelerating. 

Recently, Bali and Jain 62 have obtained a generalized expanding and shearing anisotropic 
Bianchi type I magnetofluid cosmological model for perfect fluid distribution in general 
relativity. Motivated by the situations discussed above in regard with bulk viscous cos- 
mologies, we extend their work by including an electrically neutral bulk viscous fluid as 
the source of matter in the energy- momentum tensor. The paper is organized as follows. 
In Sec. 2, we review the solutions and the main results of Bali and Jain. 62 Here we have 
shown that all the solutions obtained by Bali and Jain 62 and Bali 71 can be derived from 
our solutions. In Sec. 3, we obtain the solutions for a universe filled with a bulk viscosity 
in presence of magnetic field. Sec. 4 includes the bulk viscous cosmological solutions in 
the absence of magnetic field. In Sec. 5, we discuss our main results and summarize our 
conclusions. 

2. A Magnetic Fluid Universe- Revisited 

In this section, we review the solutions obtained by Bali and Jain. 62 We consider an 
anisotropic homogeneous Bianchi type I metric in the form given by Marder 63 

ds 2 = A 2 (dx 2 -dt 2 ) + B 2 dy 2 + C 2 dz 2 , (1) 

where the metric potentials are functions of t only. The energy momentum tensor is taken 
into the form 

T? = (p + p) Vi v j +pgl+E{, (2) 
where E\ is the electro-magnetic field given by Lichnerowicz 64 as 

E{=fi 

and 

/' /> o-;,. (4) 

Here p, p p and £ are the energy density, isotropic pressure, effective pressure, bulk viscous 
coefficient respectively and v l is the flow vector satisfying the relation 

gijvV = -1. (5) 

p is the magnetic permeability and hi the magnetic flux vector defined by 

hi = I *F jiV i, (6) 
I 1 

where is the dual electro-magnetic field tensor defined by Synge 65 to be 
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Fij is the electro-magnetic field tensor and eijki is the Levi-Civita tensor density. Here, the 
comoving coordinates are taken to be v 1 = = v 2 = v 3 and v 4 = -j. We take the incident 
magnetic field to be in the direction of ir-axis so that hi ^ 0, h 2 = = h 3 = h 4 . This leads 
to F12 = = F13 by virtue of (6). Also due to assumption of infinite conductivity of the 
fluid, we get F\± = = F24 = F34. The only non-vanishing component of F^ is F23. 
The first set of Maxwell's equation 

Fij-k + Fjk-i + F kld = 0, (8) 

leads to 

F 23 = I (const.), (9) 
where the semicolon represents a covariant differentiation. Hence 

AT 

ftl = Wc- (10) 

The Einstein's field equations with time-dependent cosmological constant 

R i - \ R 9i + A 9t = -8ttT/, (c = 1, G = 1 in gravitational unit) (11) 
for the line element (1) has been set up as 

o a2 (- I 2 \ -^44 C44 i'H 1 1 ' , .2 , 

87rA {p -2^)=— — ~ -BC + ^B + ^C- KA > (12) 



\ C44 -B4C4 

~~B C ~BC~ + AB 

,2 / - 1 ^ \ ^44 C44 A\ 2 



^ A \^WS^)'-^-^ + ^- kA - ,13) 

The suffix 4 by the symbols A, B and C denote differentiation with respect to t. Equations 
(12)-(15) are four equations in seven unknowns, A, B, C, p, p, £ and A. 
Referring to Thorn, 66 current observations of the velocity-redshift relation for extragalactic 
sources suggest that Hubble expansion of the universe is isotropic today to within ~ 30 per 
cent. 67 ' 68 Put more precisely, redshift studies place the limit 

£<0.3 

on the ratio of shear, a, to Hubble constant, H, in the neighbourhood of of our Galaxy 
today. 

Following Bali and Jain, 62 we assume that the expansion (9) in the model is proportional 
to the eigen value o\ of the shear tensor af . This condition leads to 



A = (BC) n , 



(16) 
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where n is the proportionality constant and 



1 ( 2A4 i?4 C4 

" ~ ~B ~ ~C 



ai ~ 3A V A 



I (M B4, Ci 
A I A B C 



for the metric (1). Solving these fields equations, we obtain 

A 2 



, _ sin 2 (at) 



B 2 



C 2 



N 



^2n( v^J + 1 ) 



(it \ n^fM 

sin™ (at) ( tan — 



NKl 



fed- 



■ sin™ (at) ( tan 



-L 

at \ ™vm 



where M, N and L are constants of integration and 

Wwl 2 



a 2 = 



(17) 



(18) 



(19) 



n(2n-\)M 



Hence the metric (1) reduces to the form 



, , sin (at) 

rlc^ — - i 



dx 2 



n 2 Ma 



2(1-™) 



. 2(l-n) „ 

sm ™ (at)dt 



N 



2n • ,/T7 T ''I 



sm 



"(ai) 



A' 



tan • 



at 



dy 2 



sin™ (at) 



at 

tan y 



(20) 



where 



if 2 = n(2n - 1)M. 



Using the transformations 

x = X, 
VNy = Y, 



z _ 7 
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t = T, 

the metric (20) reduces to the form 



ds 



2 _ sin 2 (aT) 
Ki 



dX 2 



K 2 ■ rp\jrp2 

(1 _ n) sin » (ai )ai 



n 2 Mo~ 



+ - 



iV 



^l 2 



1 



2^(1-/3) 



sin" (aT) 



sin" (aT) 



tan ■ 



aT 



tan ■ 



aT 



dF 2 



dZ 2 



where 



The effective pressure and energy density for the model (21) are given by 



(21) 



87rp 



K l " [gg - L 2 + (2X 2 - M) sin 2 (aT)] a 2 K? 



2(n + l) 

4sin — (aT) 



4 sin "(aT) 



A, (22) 



[K 3 cos 2 (aT)-T 2 ] a 2 ^" 

07T/9 = 2(tl+1) 5 

4sin^^(aT) 4 sin "(aT) 

where K 3 = {An + l)M. 

In the absence of magnetic field, the metric (20) reduces to 



ds 2 = K 2 T 



d X 2 - J^T^dT 2 



n 2 M~ 



+ 



I T—dY 2 



2 _ fl T~dZ 2 . 
7V2 " 



The effective pressure and energy density for the model (24) are given by 

[Jf 3 - ^ 2 + (2^2 - M)T 2 ] 1 



87175 : 



87T/9 = 



4(X 2 T 2 )^ 

(*3 - i 2 ) 



4(/f 2 T 2 )^ 



1 



- + A, 



(23) 



(24) 



(25) 



(26) 



A{K 2 T 2 )^ A{K 2 T 2 )- 

If we consider A as constant and coefficient of bulk viscosity £ as zero, we obtain the 
solutions obtained by Bali and Jain. 62 If we set A = constant, £ = and n = 1, we get 
the solutions obtained by Bali. 71 

3. Bulk Viscous Solutions in the Presence of a Magnetic Field 
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The equations with bulk viscosity can be obtained from the general relativistic field 
equations as follows. 37 To obtain these equations, we replace p by the effective pressure p 
given by (4). The expansion scalar 9 is given by 

9 = VM(n + l)K^ C0 5 flT) (27) 
sin— (aT) 

Eq. (22) can be rewritten as 

4sin^^(aT) Asm' (aT) 

For the specification of £, now we assume that the fluid obeys an equation of state of the 
form 

p = 7p, (29) 

where 7(0 < 7 < 1) is a constant. 

Thus, given £ (i) we can solve the cosmological parameters. In most of the investigations in- 
volving bulk viscosity is assumed to be a simple power function of the energy density. 34- 36 

m = Hop™, oo) 

where £ and m are constants. If m = 1, Eq. (30) may correspond to a radiative fluid. 37 
However, more realistic models 38 are based on m lying in the regime < m < |. 
3.1. Model I: (£ = &) 

When m — 0, Eq. (30) reduces to £ = £ an d hence Eqs. (28), with the use of Eq. (22), 
(27) and (29) lead to 

8.(1 + 7)P ^ **Z { n+ l)VMKF<x»(*T) + 
sin « (aT) 



A - !" [2 (K 3 - I? ) + 2K 4 M sin 2 (aT)] a 2 ^" 



4sin H ^ ii (aT) 2sin^(aT) : 



(31) 



Eliminating p(t) between Eqs. (22) and (31), we obtain 

87r£ (n + l)y/MK^ cos(aT) [K 3 - I? + (2X3 - M) sin 2 (aT) 



(1+ 7 )A : 



sin (aT) 4sin H ^ tiI (aT) 
[K 3 cos 2 (aT)-L 2 ] 7 Q 2 i^ 7 

2(ti + 1) "T" 2 , „ ) 

4sin^ i (aT) 2 sin" (aT) 



where A4 = 2n 2 — 3n — 1. 
3.2. ModeZ //; (f = £ p) 
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When m = 1, Eq. (30) reduces to £ = £ /° an d hence Eqs. (28), with the use of Eq. 
(23), (27) and (29) lead to 



871-p = 



1 



1 , _ go(»+l)y / MA: i "^ r " cos(aT) 
1 +7 (n + l) 

sin n (aT) 
2/ 



[2(X 3 -L 2 ) + 2^4Msin 2 (ar)] a 2 i^ 



2(71+1) 

4sin^^(aT) 
Eliminating p(t) between Eqs. (23) and (33), we obtain 



2 sin" (aT) 



A 



1 



-I . £ (™+l)V / M-f£" 1 2n cos(aT) 

1 + 7 T^+TT 

sin n (aT) 



(33) 



Jfi" [2(X 3 -L 2 ) + 2if 4 Msin 2 (aT)] a 2 ^™ 



2(71+1) 

4sin^^(aT) 



2 sin" (aT) 



(K 3 cos 2 (aT)- L 2 ) a 2 Kf 



2(71+1) 

Asm.— — ( a T) 



4 sin" (aT) 



(34) 



Some Physical and Geometrical Features of the Models 

We shall now give the expressions for kinematical quantities and the components of confor- 
mal curvature tensor. The scalar of expansion 9 calculated for the flow vector v l is already 
given by (27). 

The rotation u> is identically zero and the shear in the model, is given by 



2 _ [(2n- 1) 2 M cos 2 (aT) + 3T 2 ] 

~ 2(71+1) 

12sin^^(aT) 



The non-vanishing components of conformal curvature tensor are obtained as 



(35) 



C 



12 K x n [L 2 - M + 2nM + GnLVM cos(aT) + K 4 M sin 2 (aT)] 



12 



2(77+1) 

12SH1 - ^ (aT) 



(36) 



°13 — — 



^ " [L 2 - M + 2nM - GnLVM cos(aT) + K 4 M sin 2 (aT)} 



C 



2(77+1) 

12 sin"^"" (aT) 



, ^ " [T 2 — M + 2nM + K 4 M sin 2 (aT)} 



2(77+1) 

12 sin^^" (aT) 



(37) 



(38) 



The model represents an expanding, shearing and non-rotating universe in general. The 
Eq. (21) requires that n(2n —1)^0. If we choose n > 0, then p and p becomes infinite 
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at T = and T = Eq. (27) shows that 8 = oo at T = 0, 9 = at T = £ and 
= — oo at T = Thus the model expands with a big bang at T = 0, the expansion stops 
at T = and collapses at T = ^ much like the closed FRW model. The magnetic field is 
responsible for this. 

It is observed that the metric exhibits singularities at T = and T = —. At T = and 
T = 2L, the model (21) has singularity of Point type, Barrel type and Cigar type if /3 < 1, 
(3 = 1 and /? > 1 respectively as given by MacCallum 69 and Islam. 70 At the singularity 
T = 0, .922 -» oo, .933 -> oo according as ±±£ < and ^ < 0. Also at T = f , 
g 2 2 — ► oo, g 33 — » oo according as < and < 0. Thus inflationary scenario exists 
near the initial singularities T = and T = — in the model. 
The cosmic time t is given by 

t = k 1 sin- (aT)dT, (39) 



where fcj is constant. Near T = 0, sin™ (aT) ~ a~Ti i.e., t ~ k 2 Ti , where fc 2 is 
constant. Thus T — > implies i — > 0. When T — > ^ then t — > constant. Hence the model 
has a finite time span in cosmic time scale. The space time is Petrov type D when L = 
and non-degenerate Petrov type I otherwise. 

The expressions for f = ZT^ZZ™ ' § are 8 iven as follows: 

Si -a 2 sin 2 (aT) 



#i [Jf 3 cos 2 (aT) - T 2 ] - [a 2 - Asin \ (aT) ] sin* (aT) 



1 



2 cos(aT) 



(2n -l) 2 cos 2 (aT)^^ 



M 



3(n+l) 2 



(40) 



(41) 



It is observed that when T — > 0, then — > 0, which shows that the material energy is 
more dominant than magnetic near the initial singularity. 
4. Bulk Viscous Solutions in the Absence of Magnetic Field 



In the absence of the magnetic field, the pressure and density for model (24) are given 
by Eqs. (25) and (26) respectively. Eq. (25) can be rewritten as 

[K 3 -L 2 + (2K 2 - M)T 2 } 



8?r(p - £0) = 

4(K 2 T 2 )^ 

where the expansion scalar (9) is given by 

lM (n + 1) 



9 = 

{K 2 T 2 )- 

Using Eqs. (29), (30) and (43) in (42), we obtain 



1 



A(K 2 T 2 )^ 



8tt 



IP ~ ioP r ' 



n VM (n + 1) 
(K 2 T 2 )^ 



[K 3 -L 2 + (2K 2 - M)T 2 ] 
±{K 2 T 2 )^r 



(42) 



(43) 
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r - A, (44) 

A(K 2 T 2 )i 



4.1. Model I: (f = &) 



When m = 0, Eq. (30) reduces to £ = £ and hence Eq. (44) with the use of Eq. (26) 
leads to 



= 8^o(n + l)VM 2(if 3 -L 2 ) + (2/f 2 - M)T 2 
(i. 2 T 2 )^ 4(if 2 T 2 )^ 

Eliminating between Eqs. (26) and (45), we obtain 



1 



2{K 2 T 2 ) 



- (45) 



, 8< (n+l)VM K 3 - L 2 + (2K 2 - M)T 2 
(1 + 7)A = — ,„ „ ON „+i 1" - 



(K 2 T 2 )~ 
(K 3 - L 2 )j 



4(K 2 T 2 Y 



- + ■ 



4(K 2 T 2 )^ A(K 2 T 2 ) 



4.2. Model II: (f = £ p) 



(46) 



When m = 1, Eq. (30) reduces to £ = £ P and hence Eq. (44) with the help of Eq. (26) 
lead to 

8np = 



, . go(n+l)y / M 

J- t 7 ; n+i 



(AT 2 T 2 )" 

2(K 3 - L 2 ) + {2K 2 - M)T 2 



1 



4(/f 2 T 2 )^ 2(K 2 T 2 )i 
Eliminating p(t) between Eqs. (26) and (47), we obtain 

A = 



1 | Co("+l)VM 



2(^3 - L 2 ) + (2if 2 - M)T 2 



1 



4(/f 2 T 2 )^" 
(#3 - L 2 ) 



2{K 2 T 2 )- 



(47) 



(48) 



4(if 2 T 2 ) :i ^ i A{K 2 T 2 )^ 

From Eqs. (46) and (48), we observe that the cosmological constant in both models is a 
decreasing function of time and it approaches a small value as time progresses (i.e., the 
present epoch). The values of cosmological "constant" for both models are found to be 
small and positive which are supported by the results from recent supernovae observations 
(Garnavich et al.; 58 Perlmutteref al.; 55 Riess et al.; 56 Schmidt et al. 61 ) 
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Some Physical and Geometric Aspects of Models 

In absence of the magnetic field, the expansion 9 is already given by (43). The expansion 
in the model stops for large values of T. The shear in the model is given by 



2 _ (2n - 1) 2 M + 3L 2 



n _ L z - M + 2nM + 6nLVM 
^11 — „^ n+i > (50) 



M + 2nM - 6nWM 

22 



u = ± vr ■ (49) 

" + 1 , , 2(71+1) V ' 

12K 1 n (aT)^ 2 

The components of conformal curvature tensor are given by 

kf + 2nM - 
12(K 2 T 2 ) 

kf + 2nM - 
12(K 2 T 2 ) 

Thus, in the absence of magnetic field, the space-time is Petrov type D if L = and non- 
degenerate Petrov type I otherwise. For large values of T , the space-time is coformally 
flat. For n > 0, when T is large, the shear dies out. 
5. Conclusions 



^22 — ; ~ n+i j (51) 



We have obtained a new class of Bianchi type I anisotropic magnetofluid cosmological 
models with a bulk viscous fluid as the source of matter. Generally, the models are expand- 
ing, shearing and non-rotating. In all these models, we observe that they do not approach 
isotropy for large values of time t either in the presence or in the absence of magnetic field. 
The cosmological constant in all models given in Sec. 4 are decreasing function of time 
and they all approach a small value as time increases (i.e., the present epoch). The values 
of cosmological "constant" for these models are found to be small and positive which are 
supported by the results from recent supernovae observations. If we consider A as constant 
and coefficient of bulk viscosity £ as zero, we obtain the solutions obtained by Bali and 
Jain. 62 If we set A = constant, £ = and n = 1, we get the solutions obtained by Bali. 71 
Thus, with our approach, we obtain a physically relevant decay law for the cosmological 
constant unlike other investigators where adhoc laws were used to arrive at a mathematical 
expressions for the decaying vacuum energy. Thus our models are more general than those 
studied earlier. 
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